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O ' Abstract 



^ i We study the properties of black holes of mass lO'^-lO^^ GeV in models with the 

fundamental scale of gravity at the TeV. These black holes could be produced in the 
collision of a ultrahigh energy cosmic ray with a dark matter particle in our galactic 
halo or with another cosmic ray. We show that QCD bremsstrahlung and pair pro- 
duction processes are unable to thermalize the particles exiting the black hole, so a 
chromosphere is never formed during Hawking evaporation. We evaluate with HER- 
WIG the spectrum of stable 4-dim particles emitted during the Schwarzschild phase 
and find that in all cases it is peaked at energies around 0.2 GeV, with an approximate 
43% of neutrinos, 28% of photons, 16% of electrons and 13% of protons. Bulk gravitons 
are peaked at higher energies, they account for a 0.4% of the particles (16% of the total 
energy) emitted by the most massive black holes in n = 6 extra dimensions or just the 
0.02% of the particles (1.4% of the energy) emitted by a 10 TeV black hole for n = 2. 



1 Introduction 



The coexistence of the electroweak and the Planck scales has been the main motivation for 
model building in particle physics during the past 30 years. Early proposals like technicolor or 
supersymmetry could explain dynamically the hierarchy between these two scales, although 
the new physics that they suggest has been so far absent in collider experiments. More 
recent proposals offer new and very interesting possibilities. In particular, the presence of 
extra dimensions could imply a fundamental scale of gravity M/j much lower than Mp = 
G^^/^ ^ 1019 GeV iU. If Md were near the TeV region, then the electroweak scale would 
just introduce a little hierarchy problem, which could be easier to solve consistently with 
collider data [2j. 

In contrast with the usual scenario, an amusing feature in these TeV-gravity models is 
that Mf) is at accessible energies and the transplanckian regime can in principle be probed 
[3]. Actually, due to the spin 2 of the massless graviton the collision of two point-like particles 
at transplankian energies {^/s ^ Md) and large distances (6 ^ M^^) should be dominated 
by gravity. In such collision one would expect the collapse of the two particles into a black 
hole (BH) of mass M ^ ^/s with an approximate cross section 



where ru is the radius of the BH horizon. Note also that the exchange of the large momenta 
required to see string resonances [1] or other quantum gravity effects would take place here at 
shorter distances: all the details of the complete theory that describes quantum gravity will 
be screened by the BH horizon. As the collision energy increases, rii grows and the collapse 
into a BH involves larger distances, a regime where classical gravity (strongly coupled but 
with no loops [5]) should work well. 

The production of microscopic BHs at the LHC has been extensively considered in the 
literature [SI El [3, El E] . It seems difficult, however, to imagine a framework able to accom- 
modate such expectations. First of all, the mechanism that defines a consistent theory of 
gravity should also manifest below M/j while being consistent with all precision data. If, for 
example, gravity derives from string theory, the string scale Ms will be a loop factor smaller 
than Md jl]: 



with n the number of large extra dimension^. Second, even if Mr, is as low as a few 
TeV, the a/s at the parton level accessible at the LHC can not be much larger. The BHs 

^See the appendix in [TU] for the relation of Mq with 
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that one may expect there will be near the threshold M^i, where they would appear almost 
indistinguishable from a massive string mode. 

In this paper we focus on much larger mini BHs, of mass up to 10^^ GeV. They could result 
from the collision of a cosmic ray (a proton or a cosmogenic neutrino) with a dark matter 
particle in our galactic halo or with another cosmic ray. First we estimate the production 
rate of these cosmogenic BHs. Then we describe their Hawking evaporation (temperature, 
emision rate, lifetime) including in the evaluation the greybody factors for the different 
particle species. An important point that we address is the formation of a chromosphere 
or a photosphere around the BH during its evaporation. We show that, despite the large 
temperature (between 1 and 300 GeV), strong or electromagnetic interactions are unable to 
thermalize the particles exiting the BH, so their energy is not altered by this factor. Finally, 
using the code HERWIG [H] (which simulates the fragmentation of colored particles as well 
as particle decay) we derive the spectrum of stable species (protons, electrons, photons and 
neutrinos) and of bulk gravitons produced in the evaporation of a BH of given mass. Previous 
works on BH production by cosmic rays refer to the collision of a cosmogenic neutrino with 
an atmospheric nucleon [12] or a cosmic ray with a nucleon in the interestellar medium |13j . 



2 Black hole production at ultrahigh energies 

We observe a flux of cosmic rays (most of them are protons free or bound in nuclei) [2] that 
reach the earth with energies of up to 10^^ GeV. Their production and propagation induces 
a flux of (still unobserved) cosmogenic neutrinos peaked at energies near 10^ GeV [14J. For 
example, a km^ area in the upper limit of the atmosphere would receive around 10^ protons 
and 10^ neutrinoJ§ of energy between 10^ and 10^^ GeV per year. 

In addition, it is thought that 90% of the matter in our galaxy (10^^ solar masses) is 
dark: 10®^ GeV of mass in a sphere of 200 kpc, with an approximate density profil^: 

Pl.r) « , , (3) 

with i? = 20 kpc (we are at 8 kpc from the center, 1 kpc = 3 x 10^^ m). This dark 
matter would be constituted by a weakly interacting massive particle x of mass m^^ ~ 100 
GeV, although could go from 10 MeV to 10 TeV if its interaction strength goes from 
gravitational to strong [15]. 



^We take the neutrino flux in Fig. 2 of Ref. [I4] . 
^We assume a CUSP dark matter distribution [TO] . 
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Therefore, we find two types of elementary processes involving center of mass energies 
above Mr,. 



(i) The collision of a cosmic ray of energy E with a dark matter particle. Here ^/s = ^2m^E 
may go up to 10 GeV. The number of interactions to produce a BH per unit time and volume 
in terms of the dark matter density (p^ = p/m^), the flux (integrated over all solid angles), 
and the cross section is 

47r / dE a,^{s) ^ . (4) 



dtdV J ' dE 

(ii) The collision of two cosmic rays. The center of mass energy for a relative direction 6 



between the two particles is y/s = y2EiE2{l — cos 6'), which can reach 10^^ GeV. Given the 
two fluxes and the cross section, the number of interactions producing BHs is just 

-|^^16.^/d^,dE.dcos«..,Wsi„./2m|. (5) 

It is then necessary to evaluate the BH production cross section aij, where i,j label a 
proton, a neutrino or a x- We will asume that all particles except for the graviton live on a 
4 dimensional brane, that there are n fiat extra dimensions of common length, and that 
is just the higher dimensional Schwarzschild radius. 



n + 2 J \MdJ Md 

where M ^ y/s is the BH mass. It is then easy to find the cross section ai^i, (or a^y.) between 
two pointlike particles; we plot it in Fig. [T] for n = 2 and Md = 1 TeV. 

The calculation is not that simple in pu or pp collisions, as at short distances the inter- 
action involves partons. Actually, the expression for pointlike particles still holds at very 
large energies, when r// > 1 fm and the BH is large enough to contain the proton. In that 
case the neutrino interacts coherently with the whole proton and (Jp^{s) just coincides with 
a^^{s). At low s, the case extensively discussed in previous literature, the neutrino interacts 
with a single parton carrying a fraction x of the proton momentum and they collapse into a 
BH of mass M = \fxs. The cross section is 



a{xs) . (7) 



At higher energies more partons in the low x region are able to produce BHs, cTpy quickly 
grows and Eq. [7] becomes no longer reliable. A value of cTpy ~ 20 mb indicates that all 
neutrinos approaching with impact parameter smaller than the proton radius will interact 
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Figure 1: Cross sections to produce a BH for n = 2 and = 1 TeV. 

with a parton to form a BH. Obviously, increasing the energy in the colhsion we can not get 
larger cross sections unless the typical BH produced has a th similar to the proton radius 
(the regime discussed before). 

Therefore, we distinguish three regimes in the pu interaction to produce a BH. At low 
energies the neutrino interacts with a single parton, the cross section grows with s but the 
average mass M of the BH is roughly constant and close to the threshold Md. Once it 
approaches 20 mb, the cross section remains constant, but the typical BH mass increases 
with the energy. In this regime the process involves multiple scattering, in the sense that the 
BH produced in the collision will also trap spectator paitons. Finally, at higher energies (not 
in the plot) the proton becomes pointlike and the cross section grows just like a,^^. In Fig. [1] 
we have modelled a smooth transition between these regimes by discretizing the proton and 
discounting the overlapping between parton cross sections. The pp collision (also in Fig. 1) 
is completely analogous, although the 20 mb bound on app is saturated at i/i ~ 10^ GeV 
(notice that this effect is still negligible at LHC energies). 

We can now estimate the BH production rate in the two types of processes discussed 
above. In Fig. [2] we plot the number and mass distribution of BHs produced in the collision 
of a cosmic ray (a proton or a cosmogenic neutrino) with another cosmic ray or with a dark 
matter particle per year and per cubic astronomical unit (1 AU = 1.5 xlO^^ m, the mean 
earth-sun distance). We have taken the expected dark matter density at our position in the 
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Figure 2: Spectrum of BHs produced in cosmic ray collisions for n = 2, Md = 1 TeV, 
= 100 GeV, and the upper neutrino flux in [H]. 



galaxy (p = 0.3 GeV/cm^) and = 100 GeV. 



3 Black hole evaporation 

Once produced the BH will go through a fast balding phase, where it loses its gauge hair and 
asymmetries, and a relatively brief spin-down phase [8]. Most of its lifetime the BH will be 
Schwarzschild-like, and it will emit Hawking radiation [17] with an approximate black body 
spectrum of temperature 



Notice that given rn, the energy radiated by the BH will not depend (up to factors of order 
one) on the number of extra dimensions: on dimensional grounds E ~ A2+nT'^^"' ~ l/r|^ ~ 
T^. Since a bulk and a brane field see a BH of the same temperature (T is constant along 
the BH surface), the BH will emit a similar amount of energy of both fields [T8] . 

The spectrum of particles exiting the BH must cross the strong gravitational potential 
near the horizon in order to escape to infinity [19]. This effect is usually described in terms 
of the greybody factors Tg = as/A2+n, where cxs is the absortion cross section for a particle of 
spin s living in 4+n dimensions and is the BH area seen by that particle. The average 
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lifetime [s] 



Figure 3: Correlation between mass, temperature, and lifetime of a BH for Mr, = 1 TeV 
and n = 2, 6. 



emission rate for a 4-dim particle of energy uj is then 

d^ ^ A2 qF, 

being q the multiplicity of the species. Throughout the paper we use for the 4-dim particles 
the (numerical) greybody factors given in [20], together with the expressions in [21] for the 
higher dimensional graviton. The emission into bulk gravitons can be significant for a large 
number of extra dimensions, for example, it accounts for a 16% (1.6%) of the radiated energy 
for T = 10 GeV and n = Q [n = 2). In our numerical estimates we will only consider the 
emission of (relativistic) particles lighter than the BH temperature, including at T < 1 GeV 
no colored particles but pions as fundamental degrees of freedom. Notice that the emission 
will be dominated by quarks and gluons, as these cosmogenic BHs have a temperature above 

It is straightforward to integrate over all frequencies and sum over all particle species 
to deduce the BH mass loss per unit time. In Fig. [3] we plot the correlation between lifetime, 
mass and temperature for Md = 1 TeV and n = 2, 6. We see, for example, that a 10^^ GeV 
BH lives (at rest) around 10~^^ s and has an initial temperature of 0.6 GeV for Md = 1 TeV 
and n = 2. 
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4 Chromosphere around evaporating Black Holes 



An important point raised by Heckler in 1996 ^22j is that BHs above some critical temperature 
Tc ~ me/oi^^'^ may form a surrounding photospere (a plasma of electrons and photons) of 
outer temperature T ^ mg. The photosphere would thermalize to this low temperature the 
particles exiting the BH, changing dramatically the initial greybody spectrum. Moreover, 
he argued that QCD processes would also define a chromosphere in BHs of T > Aqcd, 
a fact that could affect, for example, the Page-Hawking limits on primordial BHs. After 
Heckler's initial claim there have been several analysis of photo/chromosphere formation 
[23t [23], although its existence has been considered controversial. In particular, none of the 
several codes simulating BH production at the LHC [7] has included its effect. 

Recently two different groups (Alig, Drees and Oda [25], and Carr, MacGibbon and Page 
[2S|) have reanalyzed this issue and have concluded that bremsstrahlung and pair production 
processes are not able to form a photo/chromosphere around evaporating BHs. The key 
argument is that the scattering of two particles radiated away from a BH can not be treated 
in the same way as the ordinary collision of a particle against a target, since the kinematics 
are completely different. In the radial case the particles are not coming from an infinitely 
far past, they are created in a definite point of space-time. In addition, the two particles are 
always separating (never approaching), a fact that introduces a maximum radius in which 
the process can take place. Although a general formalism to describe the radial scattering 
is not available, it is clear that the calculation of the interaction rate as 

r = {avp) (10) 

can lead to incorrect results. In this section we use the approach in [26] to show that 
the higher dimensional BHs of mass up to 10^^ GeV under analysis here do not form a 
chromosphere: quarks and gluons escape the BH (at distances around 1/Aqcd they fragment 
into hadrons, see next section) with basically the initial energy. This result is consistent with 
the detailed MonteCarlo simulation in [25] . 

To be definite we will take a very hot (LHC-like) BH, with a temperature around T = 100 
GeV, and M^i = 1 TeV for n = 3. We find that such BH emits quarks or antiquarks of 
energy E ^ 3 T with a frequency of one per r ^ 0.8/ E = O.S/'ymq (see Table 1). Once a 
quark exits the BH, it will be localized along an approximate radial distance of A ~ 1/E 
(its reduced Compton wavelength). Therefore, it can overlap significantly only with two or 
three other quarks, being the rest of quarks separated by a distance of order 1/E or larger. 

It is easy to see that the probability that the quark interacts with one of these non- 
overlapping quarks is negligible just because they are too far. The argument goes as follows. 
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n 




1/^. 


2 


16 


1.3 


3 


9 


0.8 


4 


7 


0.7 


5 


5 


0.4 


6 


4 


0.3 





175 


20 



Table 1: Average distance between consecutive electrons (l/t'e) and quarks iX/vq) in reduced 
Compton wavelength units {1/E). The last line refers to the primordial 4-dimensional BHs 
in [26] 



Suppose that the quark qi exiting the BH interacts with a quark go at a distance of A; > 1 
Compton wavelengths (1/-E) in the BH frame. The bremsstrahlung process is best under- 
stood in the rest frame of go! there qi interacts with the field generated by the static go; goes 
off-shell and emits a gluon. Since the world line of gi has a beginning, its shorter distance 
d with go will correspond to the moment when it appears in the BH horizon (except if it 
is emitted within a small solid angle of order I/7, see below). Suppose that right in that 
moment gi receives a gluon previously emited by go; it is straightforward to find that this 
gluon has been traveling a time/distancJll 

t = k 'j/E = k/rug (11) 

in the go frame. Since d ~ k/rriq ~ k/Aqcn is the minimal distance with gi that go may 
detect, the interaction will take place at typical distances where QCD is not effective (or, 
equivalently, through the exchange of momenta below the infrared cutoff Aqcd)- In [26] 
the suppression in this radial cross section with minimal distance d is estimated by the 
contribution of impact paramenters 6 > in a regular cross section (see Fig. 4). This 
suppression for the interaction between the two quarks can be also understood in the BH 
or the center of mass frames. There the time/distance that the gluon has been traveling 
is much shorter, t' ^ k/E. However, this time is already too large for the virtual gluon of 
energy of order E required by the bremsstrahlung kinematics. Notice that it is the energy 
of this virtual gluon, and not its invariant off-shellness > A^cdi "what determines the 
maximum time that it can travel without violating Heisenberg uncertainty principle. 

The causality or minimal distance argument just outhned suffices to disregard the for- 
mation of a photospere (or a chromosphere) in high temperature 4-dim BHs. The reason 



*The contributions out of the light cone are exponentiaUy suppressed in the propagator. 
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Figure 4: In a regular scattering qi comes from —oo with impact parameter 6, whereas here 
qi is created at a minimal distance d. 

is that the average distance (z/e-E")"^ (or (ugE)'^) between consecutive electrons (quarks) 
emitted by the BH is 175 (20) times larger than the typical distances dominating the QED 
(QCD) bremsstrahlung cross section. Moreover, this distance (in reduced Compton wave- 
length units) does not depend on the BH temperature, as both the number of particles and 
their typical energy grow linearly with T. In 4 + n dimensions we find (see Table 1) that 
this argument always suppresses the formation of a photosphere, but not of a chromosphere 
if n > 2. In these cases the distance between a quark and two or three other quarks is 
small (they overlap). We will then use a second argument [25], EH [27] that disfavors the 
multiple interactions required to form a chromosphere: the existence of a maximum radius, 
i^brem ~ ^/-^QCD the BH frame, where the interaction can take place. 

This maximum radius appears in the interaction of two particles moving with a relative 
angle larger than I/7; if they move in the same direction their distance may not increase, 
but the density in their center of mass frame is diluted by a Lorentz factor and becomes 
too low to give interactions [26] (6' < I/7 defines an exclusion cone). The key observation is 
that each quark can complete at most one bremsstrahlung interaction before crossing Threm, 
so after that interaction its (tranverse) distance with the particles out of the exclusion cone 
will be much larger than 1/Aqcd- To understand that, let us suppose that, right when it is 
created next to the BH horizon, a quark q absorbs a virtual gluon of ~ and goes off 
shell. In the q rest frame both its energy and its off-shellness after absorbing the gluon are 
of order A, which in a QCD bremsstrahlung should be just above Aqcd- Then the virtual q 
lives a time of order 1/A and decays into the final quark and gluon. Now, going back to the 
BH frame we observe the lifetime of the virtual gluon a Lorentz factor larger, so the typical 
distance that it travels will be of order 7/A ^ 1/Aqcd ~ "'^t-rem- Except for small values of 
7 {i.e., non-relativistic quarks emitted by BHs of temperature close to Aqcd) this argument 
establishes that quarks exiting the BH cannot interact with each other more than once, as 
would be necessary to form a chromosphere. Our conclusion agrees with the simulation of 
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Figure 5: Spectrum of stable particles and bulk gravitons (dashed) at T = 10 GeV for 
Md = 1 TeV, n = 2 (left) and n = 6 (right). 



BH production and evaporation at the LHC in [25] . 



5 Spectrum of stable particles 

The cosmogenic BHs under study here are produced at astrophysical distances from the 
Earth, therefore, unstable particles resulting from their evaporation have plenty of time to 
decay. A BH emits stable neutrinos, electrons, photons and gravitons, but mostly, it emits 
quarks and gluons that will fragment into hadrons and then shower into stable particles. In 
this section we evaluate the total spectrum of stable species from a BH of mass up to 10^^ 
GeV. Our results are analogous to the ones obtained by MacGibbon and Webber in [28] 
for primordial BHs [29j. Of course, we use an updated MonteCarlo jet code (HERWIG 6.5 
|llj ) and include the effects of the extra dimensions, namely, appropriate greybody factors 
and bulk graviton emission. In addition, while in primordial BHs the spectrum corresponds 
to a given temperature (T changes only within astrophysical time scales), here we evaluate 
the total spectrum resulting from the complete evaporation of the BH, which includes a 
(relatively small) contribution from the high temperatures briefly reached at the end of its 
lifetime. 

Several comments are here in order. As we have explained in the previous section, when 
the quarks and gluons exit the BH their probability of interaction is small. We will then 
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Figure 6: Spectrum of stable particles and bulk gravitons (dashed) from a BH of M = 10 
TeV (left) and M = 10^° GeV for Mn = 1 TeV and n = 2. 

assume that the jets that they define are similar to the ones produced in e+e~ collisions to 
qq or (the fictitious process) gg. We have used HERWIG [H] to simulate the jets produced 
by any particle that is light at a given temperature, including massive gauge bosons and the 
top quark (but not the Higgs boson nor a dark matter particle). Finally, we have added 
together the number of particles and antiparticles (they are produced at the same rate) and 
the three (Majorana) neutrino species (at astrophysical distances their flavour oscillates). 

In Fig. 5 we plot the power spectrum emitted by a BH of T = 10 GeV for n = 2, 6. 
We observe that it is dominated by energies around Aqcd, although the primary greybody 
spectrum peaks at 30 GeV. This is manifest in the flux of gravitons, since they decouple 
(their number is not increased by the showering of unstable particles). At high energies 
it is possible to distinguish the primary greybody spectra for some of the particle species. 
At this temperature the particles emitted onto the brane consist of an approximate 43% of 
neutrinos, a 28% of photons, a 16% of electrons and a 13% of protons. The bulk gravitons 
are a 0.02% of the total emitted particles for n = 2 or a 0.4% for n = 6. Their typical energy 
is higher, so they account for a 1.4% (n = 2) or a 16% (n = 6) of the total energy radiated. 

In Fig. 6 we give the total spectrum from BHs of mass 10 TeV and 10^*^ GeV for n = 2. 
These masses correspond to initial temperatures of 120 and 1.2 GeV, respectively. We find 
that the spectra are dominated by the emission at these initial temperatures, in particular, 
90% of the energy is emitted by the 10^° GeV BH when its temperature has only increased 
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from 1.2 to 2.6 GeV (see Fig. [3]). The dominance in the spectrum of energies around 0.2 
GeV and the approximate fraction of different species are features that depend very mildly 
on the BH mass, as it is apparent in these plots. Of course, in smaller BHs of higher initial 
temperature the relative weight of high energies in the spectrum is larger. 

6 Summary and discussion 

Models with the scale of gravity at the TeV must confront the fact that in nature there are 
processes at much higher energy. In particular, in the collision of two cosmic rays the center 
of mass energy can go up to 10^^ GeV. At these energies for impact parameters short enough 
the particles will be trapped inside a gravitational horizon and form a mini BH. We have 
estimated the production rate of these cosmogenic BHs (in Fig. 2) and have analyzed their 
properties. 

(i) We have found that their lifetime goes from 10~^^ s for a light BH of M = 10 TeV to 
10~^^ s for M = 10^^ GeV. In Fig. |3]we plot the correlation between BH mass, temperature 
and lifetime for n = 2, 6. 

(ii) We have shown that bremsstrahlung interactions between the particles exiting the BH 
are unable to form a chromosphere. Although the average separation between these particles 
(in Table 1) can be in some cases small, a QCD process requires a typical time to develop, 
and this time is such that the particle can not interact more than once at r < 1/Aqcd in 
the BH frame. 

(in) The greybody spectrum emitted by the BH onto the brane is dominated by colored 
particles. These form jets that result into a spectrum of stable particles peaked at Aqcd- 
We obtain an approximate 43% of neutrinos, a 28% of photons, a 16% of electrons and a 
13% of protons. These two features in the spectrum depend very mildly on the number of 
extra dimensions (Fig. 5) or the BH mass (Fig. 6). 

(iv) Bulk graviton emission is relatively larger for high BH masses and large number of 
extra dimensions, and it is peaked at higher energies than for the rest of species. It accounts 
for a 0.4% of the total particles (16% of the energy) emitted by a 10^° GeV BH for n = 6 or 
just a 0.02% (1.4% of the energy) for M = 10 TeV and n = 2. 

We find remarkable that, even though there is no chromosphere, the spectrum of stable 
particles from the evaporation is equally peaked at low (0.2 GeV) energies: the spectrum 
provided by one of these BHs at astrophysical distances would not be too different whether 
there is or there is not a chromosphere around it, in both cases the scale is fixed by Aqqd- 
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Notice, however, that the signal at the LHC in one case or the other would be clearly different 
(10 jets of 100 GeV each versus an expanding shell of 1000 hadrons). 

We think that the analysis presented here is a necessary first step in the search for possible 
observable effects from BH production by cosmic rays. 
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